Electric-field-driven phase transition in vanadium dioxide 
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We report on local probe measurements of current-voltage and electrostatic force-voltage char- 
acteristics of electric field induced insulator to metal transition in VO2 thin film. In conduct- 
ing AFM mode, switching from the insulating to metallic state occurs for electric field threshold 
£ ~ 6.5 x 10 7 Vm _1 at 300 K. Upon lifting the tip above the sample surface, we find that the 
transition can also be observed through a change in electrostatic force and in tunneling current. In 
this non-contact regime, the transition is characterized by random telegraphic noise. These results 
show that electric field alone is sufficient to induce the transition, however, the electronic current 
provides a positive feedback effect that amplifies the phenomena. 

PACS numbers: 72.70.+m,72.20.-i,71.30.+h 



h 



I 

C 

O 

o 



> 
in 



o 



X 



Understanding the non-equilibrium transport prop- 
erties of strongly correlated systems such as Mott 
insulators |T] is a fundamental issue of condensed mat- 
ter physics. This out of equilibrium many body prob- 
lem presents a formidable challenge to theory further- 
more, the non-linear current voltage characteristics I(V) 
observed in correlated materials are of interest for the de- 
velopment of new electronic components; as exemplified 
recently with the discovery of the memristorsQ. 

Remarkable non-equilibrium transport properties due 
to electric field induced switching from insulating to a 
conducting regime have been observed in manganitesQ, 
vanadium dioxide (V02)@41] and magnetite thin films 
or nanoparticles[9|. 

The material VO2 is particularly interesting as it un- 
dergoes a first-order thermal induced Insulator to Metal 
Transition (IMT) at 340 K^U^ where the resistivity 
changes by several orders of magnitude. The insulat- 
ing state can be viewed as a solid of singlet vanadium 
atoms dimers in the Heitler-London limit, implying elec- 
tron correlations llj . The electrical field induced break- 
down of this insulating state led to large current jumps 
(switching) in the I(V) curves^ |(|. 

The most important issue about this non-equilibrium 
transition is to understand the role of the electronic cur- 
rent. In non-equilibrium transport models based on elec- 
tric field assisted tunneling such as in Poole-Frenkel or 
Zener mechanisms, only the electric field matters and the 
electronic current plays no role. In these models, the am- 
plitude of the current can only increase exponentially fast 
with applied voltage. The switching phenomena - which 
led to a faster than exponential increase of the current - 
is puzzling as it implies either a many-body effect, such as 
a phase transition between an insulating and a metallic 
state, or some other non-linear effects due to avalanche 



breakdown or Joule heating. While some experiments 
suggest that Joule heating could be responsible for the 
transition 0, others suggest this is negligeable[l2j|. 

To answer that question, we used an Atomic Force 
Microscope (AFM) to perform local measurements on a 
VO2 film in the conducting AFM mode and in the Elec- 
trostatic Force Microscopy (EFM) mode. In the conduct- 
ing AFM mode, the tip is in close contact with the sample 
and switching in the I(V) curve is observed. When the tip 
is retracted, we find that switching events can also be ob- 
served in the electrostatic force - voltage F(V) character- 
istics and in the tunneling current-voltage characteristics. 
In these last two modes, where no large current flows, the 
data present Random Telegraphic Noise (RTN) , indicat- 
ing that the system is bi-stable, jumping randomly from 
the insulating to the conducting state when an electrical 
field is applied. This shows that electric-field induced 
breakdown of a correlated insulating state is possible 
without any electronic current. However, the observa- 
tion of RTN in those data, not observed when a large 
electronic current flows across the sample, suggests that 
the current plays a role in the transition through a posi- 
tive feedback effect. 

A 80 nm VO2 thick film was deposited on a Si doped 
substrate by pulsed laser ablation [13[ . AFM microscopy 
gives RMS surface roughness about 5 nm as shown in 
the insert of Fig. QJ). Figure QJi shows the sample resis- 
tance as the temperature is swept across the IMT. Local 
probe measurements are carried out with a home-built 
variable temperature AFM installed in Ultra High Vac- 
uum. The force detection is based on measurement of 



frequency shift of a tuning fork 14]. Typical amplitude 
oscillation of the tuning fork is about 5 nm. The Pt/Ir 
tips glued onto the tuning fork were prepared by electri- 
cal etching or by mechanical cutting. 
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Figure Q] shows typical I(V) curves measured with the 
tip in contact with the surface, at different tempera- 
tures. At 300 K, the first current jump is observed at 
voltage bias of 5.2 V, which corresponds to an electric 
field threshold £ = 6.5 x 10 7 VmT 1 . As the tempera- 
ture is increased, this threshold is shifted to lower volt- 
age; at T = 330 K, the first switching voltage is 4.5 V, 
i.e. £ = 5.6 x 10 7 VmT 1 . These electric field thresh- 
olds are in good agreement with previous local probe 
measurements 0, HJ of switching current in VO2. In con- 
trast to measurements performed with macroscopic elec- 
trodes, our data can present several switches per I(V) 
curve. Probably, each additional current jump in the 
I(V) curve corresponds to the spreading of the transition 
area below the tip. 

For temperature T = 360 K, no current jumps are 
observed as expected for the sample in the metallic state. 
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FIG. 1. (color online) Panel (a) Temperature dependence of 
the in-plane resistance of the VO2 thin film deposed on doped 
silicon; Panel (b) I(V) curves measured at 300 K, 330 K, and 
I(V)/10 measured at 360 K. Insert, AFM topography of the 
sample. Panel (c) Zoom on the current jumps, indicated as a 
dashed box on panel (b). 



With this setup, one can now lift the tip and perform 
electrostatic force measurements. If the tip is far enough 
(more then 25nm), no current can be detected; however, 
the electrostatic force causes a shift of the tuning fork 
resonance frequency. Given that the electrostatic energy 
of the tip-sample capacitor C is U e i — (I/2C)A</> 2 , where 
A(f> = (ptip — (('sample + V is the electrochemical potential 
difference between the tip and the sample added to the 
applied voltage V, the electrostatic force is IBl 16j|: 



F 



ldC 
2~dz~ 




-4 4 8 

Voltage (V) 

FIG. 2. Panel (a): Tuning fork frequency shift A/ vs. tip- 
sample voltage difference V for the tip located at 35nm above 
the surface. It shows the characteristic quadratic dependence 
of the electrostatic force with voltage. Small steps appear at 
high voltage, ~ 7 V, as shown in the inset with the parabolic 
background removed. They are the consequence of electric 
field induced IMT. 



with the effective tip-sample distance z = d + Xd, 
where d is the tip-sample distance and An the Debye 
Screening Length (DSL). Following Ref. [17|, the fre- 
quency shift A/ is : 



dF 1 d 2 C . l2 
A/ = -^=2^ 



(2) 



The experimental quadratic behavior of the frequency 
shift vs. voltage is shown in Figure [2] 

At tip-sample distance d = 35 nm and high voltage 
~7y, one observes small jumps in the electrostatic force 
A/ vs. V curve. According to formula [T] these jumps 
could be interpreted either as changes in chemical po- 
tential or as changes in tip-sample capacitance. As both 
the chemical potential and the DSL should be modified 
at the transition from the insulating to the conducting 
state, an electric field induced IMT provides a natural 
explanation for the observed jumps [18]. As a change in 
the DSL induces a change in the effective tip-sample dis- 
tance z, its effect on the frequency shift can be estimated 
from SAf/Af — —Sz/z obtained from the derivation 
of Eq[TJ The magnitude of the jumps are in the range 
0.01 — 0.05 Hz. As the full frequency shift due to elec- 
trostatic force is ~ 5 Hz, the relative change is about 
0.2 — 1 %, which implies, with the tip-sample distance 
d = 35 nm, a change of the DSL about ~ 0.07 — 0.35 nm. 

The DSL for VO2 can be estimated from the carrier 
density n : 
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At the thermal-tuned IMT, the electron density changes 
by several orders of magnitude. In the metallic state, 



3 



N 



200 



X 



(a) 



n 1 r - 



n 1 1 — 

Pt300K 



-VO, T>T 370K 



-VO, T<T , 300K trace 1 



w """" V0 2 T<T U| 300K trace 2 

_J I 1 I I I I I i l_ 



400 800 1200 1600 2000 

time (s) 




-40 -30 -20 -10 10 20 30 40 
AZ'(mHz) 

FIG. 3. (color online) Panel (a): Tuning fork frequency shift 
vs. time. From bottom to top, the first two curves are taken 
at two different places on the VO2 sample at 300 K; the third 
curve shows a time trace at 370 K in the metallic state, and 
finally, the upper curve is a time trace on a simple Pt thin 
film at 300 K. Panel (b) Histograms of all time traces (using 
the same color code as panel (a)). Curves are normalized to 
set their maximum to 100 counts for each histogram. Insert 
panel (b): two-states fiuctuator model discussed in text. 



Hall conductance measurements [l9j give a carrier density 
about n = 10 22 cm -3 , while capacitance measurements 



give n — 10 19 cm a 20]. Using the dielectric coefficient 
e = 30 [20], this led to a DSL in the range \ D ~ 0.07 - 
2 nm. In the insulating state, capacitance measurements 
give a carrier density ranging from n = 5 x 10 15 cm -3 to 



5 x 10 18 



20] . This leads to a DSL in the range 



Xd — 3 — 90 nm. With those values, we find that the 
change in DSL at the thermal tuned IMT is in the range 
[1 — 90 nm]. Thus, a change of DSL of one nanometer 
due to the electric field induced IMT is quite possible and 
provides a plausible explanation for the observed jumps 
in the frequency shift. 

In contrast to the contact mode, we find that the 
electrostatic force detected transition is characterized by 
non-Gaussian noise in the frequency shift measured as 
function of time. Many time traces, for fixed tip posi- 
tion 35 nm above the surface, lasting about one hour, 
have been taken which all present RTN [21] as shown 
by the bottom curve of Fig. [3k,. A histogram of each 
time trace, Fig. [3p, shows that this noise is obviously 
non-Gaussian at 300K with two peaks corresponding to 
two stable states. We checked that the observed RTN 
could not arise from the instrument. First, RTN is never 
observed when measurements are performed on simple 



metal Pt surface, see Fig. [3^, and second, RTN disap- 
pears when the temperature is increased at 370K above 
the IMT transition of VO2. In those last two cases, 
we find instead that the noise has a Gaussian distribu- 
tion, Fig. RTN such as observed here is the charac- 
teristic signature of a two-states fiuctuator. Two-states 
fluctuators can emerge from localized states of nanoscale 
devices or nanoparticles|22j, where fluctuations between 
the two states is due to trapping-detrapping of single 
charges. However, two-states fluctuators can also emerge 
from systems with many degrees of freedom. Typical 



examples are provided by Josephson junctions 23] or su- 
perconducting nanowiresf^jj] where the system switches 
from the superconducting to the normal state when the 
superfluid current exceeds some critical value. Further- 
more, RTN has also been observed in transport mea- 
surements of strongly correlated material such as un- 
derdoped cuprates[2£ 



and manganitcs 26], interpreted as 
two-levels fluctuators emerging in systems with many de- 
grees of freedom. In our system, the two-states fiuctuator 
is built up from the two macroscopic states of VO2 : one 
is insulating, one is conducting. 

It has a double potential well structure 



2G. 27, as 



drawn schematically in the insert Fig. [3Jd. RTN is pro- 
duced as thermal excitations drive the system randomly 
in one state or the other. In this model, the average time 
spent in one state i is given by the Arrhenius law[28j : Tj = 
To exp -j^kp where E^ is the energy of the state i and r is a 
microscopic constant which depends on the thermal cou- 
pling to the external bath. Using this formula, one finds 
the energy difference Ejj — En = fc^Tln^ = 17 meV, 
where the occupation times tjj = 5.5 s and td = 2.9 s 
are obtained [2^] from counting the occupation times in 
the lowest time trace Fig. [3^. 

Note that this energy level difference should depend 
on temperature and applied electric field. Furthermore, 
the detailed characteristics of the noise may change 
from one sample location to another one, as exempli- 
fied by the time trace 1 taken on another location of 
the sample, shown Fig. [3^. Interestingly, recent theoreti- 
cal works [29j have shown that electronic inhomogeneities 
such as stripes could lead to RTN, whose detailed char- 
acteristics could reveal whether such noise is due to cor- 
related fluctuations rather than independent switcher. 
However, such an analysis will require a more detailed 
spatial mapping of RTN than presented here. 

Figure [4] shows that RTN only appears above some 
electric field threshold. For a tip-sample distance of 
35nm, RTN disappears for voltage values lower than 
~ 4 V. On the other hand, for a fixed voltage value 
~ 8 V, RTN disappears for tip-sample distance larger 
than 90 nm. Using for the DSL the value estimated 
above, we find that the electric field threshold is of the or- 
der of £ ~ 10 s Vm^ 1 , which is comparable to the electric 
field threshold found for current switching in the contact 
mode. 
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FIG. 4. (color online) EFM time traces for different electrical 
field values at 300K. Panel (a), the tip/sample distance is set 
to 35 nm and voltage is changed from 8V to 0V. Panel (c), the 
tip/sample voltage is set to 8 V and the tip/sample distance 
is varied from 30 nm to 130 nm. Panels (b) and (d) show the 
histograms of each time traces presented in panels (a) and (c) 
respectively. Non-Gaussian noise is only found for tip/sample 
distance below 90nm and tip/sample voltage above 4V. 



of the current in the transition. When the electronic cur- 
rent is allowed to flow, it induces some positive feedback 
- possibly through Joule heating - and drives the system 
deeper into the metallic state. On the other hand, if the 
current is not allowed to flow, as there is only a small free 
energy difference between the insulating and conducting 
states, the system keeps switching between the two states 
on a large electric field range. 

To conclude, we have presented conducting AFM and 
EFM on VO2 thin him. In the contact mode, current 
switching is seen in the I(V) curves. With no Joule heat- 
ing (tunnel mode and EFM mode), the system does not 
have any positive feedback and RTN is observed in the 
data. The observation of these fast switching events in 
the absence of Joule heating indicates that many-body 
effects should be at the origin of this phenomena. Such 
a many-body phenomena could involve an electric field 
induced IMT or possibly some many-body versions of 
Zener tunneling as recently studied theoretically 0. 

A.Z. and H.A acknowledge support from ANR grant 
09-BLAN-0388-01. R.G, Y.L. and R.L. acknowledge sup- 
port from UNC-Chapel Hill Institute for the Environ- 
ment Carolina Energy Fellows Program. 
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FIG. 5. (color online) I(V) curve in the tunnel regime. 



An intermediate situation between good tip-sample 
contact, shown Fig. [TJ and purely electrostatic mode, 
shown Fig. [3J is the tunnel regime. In this regime, RTN 
is also observed in the current-voltage characteristics, 
shown Fig. [SJ At low voltage, the sample is in the in- 
sulating state and at high voltage, the sample is in the 
conducting state. However, on a wide intermediate elec- 
tric field range, the system switches randomly from the 
metallic to the insulating state, producing the observed 
RTN. 

The observation of RTN in both the electrostatic force 
and tunnel current but not in the I(V) curves measured in 
the contact mode point toward an explanation of the role 
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